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I. INTRODUCTION 

K13 decays, along with the leptonic decay of the kaon, 
are the gold-plated channels for a precise determina- 
tion of | V us |, where V us is the element of the Cabibbo- 
Kobayashi-Maskawa matrix (for recent reviews see 
Q). The amplitude of the process involves the matrix 
element of the strangeness-changing vector current be- 
tween a kaon and a pion, written as: 



(n°(p')\s^u\K+(p)) 
1 



V2 



[( p '+pW + (*) + (p-pV-(*)]. (i) 



where /+ (t) is the vector form factor and the combination 



/o(t) = /+(*) + 



Ml 



/-(*) 



(2) 



is known as the scalar form factor. The matrix element 
for the charged pion and the neutral kaon is related to 
([2]) by isospin symmetry. 

The K13 decay rates were measured for the four modes 
(K = K ± l K° and I = fi, e) by several experimental 
groups [f|-[lH. The rates are conveniently written as 



Ki 



G%M 5 K 2 



192tt 3 



K°EW 



f + (0)V us I e K (1 + A) , (3) 



where Gf is the Fermi constant, Ck is the Clebsh- 
Gordan coefficient equal to 1 (I/a/2) for the neutral 
(charged) kaon decays, 5ew is a short-distance elec- 
troweak correction, and A accounts for the electromag- 
netic and isospin-breaking corrections. The form fac- 
tors enter through the value /+(0) at zero momentum 
transfer and the phase space integral I l K , which depends 
on the shape of the form factors in the physical range 
Mf < t < (M K - NU) 2 . 

For a precise determination of \V US \, it is important 
to improve the accuracy of the parameterizations of the 



form factors using additional theoretical and experimen- 
tal information. Thus, the low-energy theorems based 
on chiral symmetry provide values of the form factors at 
some special points inside the analyticity domain |16j |- 
[2lj . On the unitarity cut, which extends from t+ = 
(M K + M T ) 2 to oo, the Fermi- Watson theorem [H, HU 
implies that, below the inelastic threshold t- m , the phase 
is available from the corresponding partial wave of pion- 
kaon elastic scattering. Furthermore, recent measure- 
ments of t — V Ktxv t decays [HI provide experimental 
information also on the modulus in the same region. 

Analyticity is the ideal tool for relating the information 
from the unitarity cut to the semileptonic range. Sev- 
eral comprehensive dispersive analyses were performed 
recently, using either coupled channels Muskhelishvili - 
Omnes equation s [25| -[28| or a single-channel Omnes rep- 
resentation [29l . l30l |. The dispersive representations can 
be extrapolated below the cut, providing information on 
the shape of the form factors in the region. How- 
ever, direct applications of the dispersion relation for the 
data analysis are not usual , although exceptions are the 
Omnes- type relations 0, [3(| , used recently in the data 
analyses by NA48 0, KLOE [ij] and KTeV [3 Collab- 
orations. Such an analysis of BELLE data for r decays 
that probe the vector form factor is found in [3lj . 

The purpose of the present paper is to discuss the im- 
plications of analyticity for the phenomenological anal- 
yses using an alternative approach proposed some years 
ago [32H34j |. known as the method of unitarity bounds. 
We use the fact that a bound on an integral involving the 
modulus squared of the form factors along the unitarity 
cut is known from the dispersion relation satisfied by a 
certain QCD correlator. Standard mathematical tech- 
niques then allow one to correlate the values of the form 
factor or its derivative at different points. For the Kir 
form factors, the method was applied in [35U4ll ]. The 
latest applications 0, |4l| led, in particular, to stringent 
constraints on the shape of the scalar form factor at low 
energies. 

In the present work, we consider both the vector and 
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the scalar form factors and focus on the phenomenologi- 
cal consequences of analyticity and unitarity for anal- 
yses. One of our aims is to present simple analytic con- 
straints, easily implementable in phenomenological stud- 
ies, on the free coefficients of the parameterizations used 
in fitting the data. In contrast to other recent works, 
the present work does not require any input about the 
absence of zeros on the real energy line or in the complex 
energy plane, nor does it require any knowledge of the 
phase of the form factor in experimentally inaccessible 
regions. Thus, the results of the present work are a rig- 
orous consequence of the general principles and do not 
have any model dependence. 

We start by giving in Sec. [TT| a brief review of our 
theoretical framework. In Sec. Mil we present in de- 
tail the input quantities used in the application of the 
method to the K13 form factors. In Sec. IIVI we review 
the main parameterizations used in the analyses, 
emphasizing the merits and the shortcomings of each of 
them. In Sec. [V] we concentrate on the parameteriza- 
tion based on the standard Taylor expansion at t = 
and derive explicit constraints on the coefficients of the 
expansion. To facilitate further applications, we present 
the results as simple quadratic expressions for arbitrary 
input values of the form factors at special points inside 
the analyticity domain (the origin t = and the Callan- 
Treiman(CT) point). For numerical illustrations, we use 
as input the precise values obtained recently from cal- 
culations in ChPT and on the lattice. We work in the 
isospin limit, but briefly discuss also the effects of sym- 
metry breaking in Sec. IV Al Further, in Sec. I VII we 
investigate the truncation error related to the higher or- 
der terms in the expansion, and in Sec. I VIII we show that 
the method allows one to derive in a rigorous way the 
domains in the complex i-plane where zeros of the form 
factors are excluded. Sec. IVIIII contains some final re- 
marks and our conclusions. 



II. FORMALISM 

A review of the formalism of unitarity bounds-in the 
standard version and the modified forms that include ad- 
ditional information on the unitarity cut- was given re- 
cently in [4l|. Here we shall present for completeness 
the approach proposed in [42| , which will be used in the 
applications made below. 

We shall denote generically the form factors /+ (t) and 
fo{t) by a function F(t), which is real analytic in the 
complex i-plane except for the unitarity cut along the 
positive real axis from the lowest unitarity branch point 
t + = (Mk + M^) 2 to 00. The dispersion relation satis- 
fied by the correlator of two strangeness-changing vector 
currents (see details in Sec. IIIip implies an inequality of 
the type- 



dt p{t)\F{t)\ 2 <I, 



(4) 



where the weight function p(t) > and the quantity / 
are known. 

According to the Fermi- Watson theorem [22|, [H| , be- 
low the inelastic threshold fi n the phase of F(t) is equal 
(modulo 7r) to the phase S(t) of a partial wave of ttK 
elastic scattering. Thus, we can write 

F{t + ^e) = \F(t)\e im , t+ < t < t in , (5) 

where S(t) is known. We define the Omnes function 



0(t) = exp 



6(f) 



dt'- , 



(6) 



where S(t) is known for t < i; n -and is an arbitrary func- 
tion, sufficiently smooth (i.e., Lipschitz continuous) for 
t > tin- From ([S]) and ([5]) it follows that the function 



h(t)=F(t)[0(t)}-\ 



(J) 



has a larger analyticity domain, namelythe complex t- 
plane cut only for t > t- m . We further assume that a reli- 
able parameterization of the modulus \F(t)\ is available 
on the same range t + < t < t; n , such that the quantity 



/' = /-/ dtp(t)\F(t)\ 2 (8) 
Jt+ 

is known. Then from Q we obtain an L 2 norm condition 
' dtp(t)\0(t)\ 2 \h(t)\ 2 < I' (9) 

tin 

for the function h(t) analytic in the t-plane cut for t > t\ n . 
As shown in [42j , ([9j can be brought into a canonical form 
by making the conformal transformation 



i(t) = 



(10) 



which maps the complex i-plane cut for t > i; n onto the 
unit disk in the z-plane defined by z = z(t). Then, Q 
can be written as 



1 

2^ 



2tt 



d%(exp( ? 0))| 2 </', 



where the function g(z) is defined as 

g(z) = w(z)u J (z)F(t(z))[0(z)}-\ 



(11) 



(12) 



In this relation, t(z) is the inverse of z = z(t), for z(t) 
defined in (|10p and w(z) is an outer function, i.e., a func- 
tion analytic and without zeros in \z\ < 1, whose modu- 
lus on the boundary is related to the weight p(t) and the 
Jacobian of the transformation (flT)]) by 



|w(exp(i6»))| 5 
2^ 



= p(t(exp(i0))) 



d*(exp(i0)) 
d9 



(13) 
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In general, an outer function is obtained from its modulus 
on the boundary by the integral 



w(z) = exp 



27T ./0 



d0- 



+ z 



lnHe ie )| 



(14) 



but in simple cases one can obtain an analytic form (see 
Sec. MI Dp . Further, the function O(z) is defined as 



O(z) = 0(t(z)), 



(15) 



and 

u(z) 



exp 



\Ain - t(z) 



dt'- 



ln|C(f) 



tin 



Vt^Uni? - t(z)) 



(16) 

From the definition (fl"2j). taking into account (JT]), it 
follows that g{z) is analytic within the unit disc \z\ < 1. 
The relation (fTTj) leads to constraints on the values of 
g and its derivatives at various points (mathematically, 
this is known as the Meiman problem). In the general 
case, consider the first K derivatives of g(z) at z = and 
the values at other N interior points z n : 



1 d k g(z) 



k\ dz k 



z=0 

sO™) 



9k, 0<k<K-l; 

z n ^0, l<n< N, (17) 



where g^ and are given numbers. Then the following 
determinantal inequality holds: 



6 



6 



,2K 



K 



(ziZ 2 



3c 



6 

[z\z 2 ) 

1 - ZlZ 2 

(z 2 ) 2K 



1 - ZiZ 2 



where 



(ziZjv) (z2Z/v) 



K 



1 - ZiZat 1 - Z 2 Zn 



(z iZn ) k 

1 - ZlZ N 

(z 2 z N ) K 
1 - z 2 z N 



,.2K 
Z N 

i-4 



>0, (18) 



/=/' 



K-l 

k=a 



(19) 



All the principal minors of the above matrix should also 
be nonnegative [H, [35[ . 

The entries of the determinant (|18p are related, by 
(JUJ), to the derivatives F^(0), j < K - 1 of F(t) at 
t = 0, and the values F(t(z n )), respectively. It should 
be noted that (fT5)l covers also the case of values given at 
complex points: if one of the numbers t n is complex, also 
the complex conjugate appears, say t n+ i = t*, and we 
have F(t(z n+ i j) = F*(t(z n j) due to the reality property. 
The same relations hold for the corresponding points in 
the z-plane and the values g(z n ). This ensures the reality 
of the determinant appearing in the inequality (TH 



We note that the formalism presented above exploits in 
an optimal way the relation (O , which is a consequence of 
the inequality (111) and of the relations ©-([SI). The stan- 
dard approach [33. [33ll35l. l37l. 38], which does not include 
additional information on the form factors on the cut, 
is obtained formally from the above relations by setting 
tin t+i when both the Omnes function 0(t) and the 
function lo(z) become unity. The implementation of the 
phase condition §5§ on the elastic cut, together with the 
relation (jl]) , involves the solution of an integral equation 
of Fredholm type (for details and references, see |4l[). 

We stress that, while © is a necessary condition fol- 
lowing from the original inequality Q , it is not sufficient 
for the fulfillment of ((4]), i.e., functions that satisfy ((9]) 
and do not satisfy Q may exist. In principle, both condi- 
tions must be imposed in order to restrict the allowed do- 
main of the parameters of interest: each condition leads 
to an allowed domain, the final region being the intersec- 
tion of the corresponding domains. As shown in [4l| . by 
a conservative choice of ij n , the results obtained from © 
satisfy also the condition ((4]). We shall place ourselves in 
this framework here. Finally, we recall that in the def- 
inition (JBJ) of the Omnes function, the phase for t > t m 
is not specified and can be parameterized in an arbitrary 
way. As shown in [4l| , the results are independent of the 
form chosen, in particular of the phase at asymptotic en- 
ergies, provided that the parametrization is sufficiently 
smooth (i.e. Lipschitz continuous). 



III. APPLICATION TO THE A i3 FORM 
FACTORS 



In this section we apply the above formalism to the 
form factors f+(t) and fo(t). We first write down disper- 
sion relations for suitable QCD correlators and show how 
they lead to an inequality of the type ((U). Then we briefly 
discuss the low-energy theorems and the information on 
the phase and modulus on the elastic part of the cut used 
as input. For completeness we give also the explicit form 
of the outer functions defined above. We work in the 
isospin limit, adopting the convention that Mk and 
are the masses of the charged mesons. A few comments 
about isospin breaking effects will be made in Sec. IV Al 



A. QCD correlators 

We consider the correlator of the strangeness-changing 
hadronic current V v = s^fu: 

W v {q) = i Jd 4 xe lqx (0\T{V^(x)V ,/ (0y}\0) (20) 

= (-s"V + 9V)ni(g 2 ) + q^n ( q 2 ). 

In QCD, the invariant amplitudes satisfy subtracted dis- 
persion relations. More exactly, it is convenient to define 



4 



the functions [H, HE [H| 1 

1 8 2 

xAQ 2 ) 



2d{Q 2 ) 2 



[g 2 n!(-Q 2 )] , 



xo(Q 2 



[Q 2 n (-Q 2 )] 



which satisfy the dispersion relations 



xi(Q 2 ) 



l 



dt 



flmlli(i) 



Wo (i + Q 2 ) 3: 



(21) 



(22) 



(23) 



and 



Xi(Q 2 ), Io=Xa(Q 2 ). 



(30) 



We evaluated these expressions taking Q 2 = 4 GeV as 
in [H [al , m s (2 GeV) = 98 ± 10 MeV, m u (2 GeV) = 3 ± 
lMeV Q and a s (2GeV) = 0.308 ± 0.014, which results 
from the recent average a s (m T ) = 0.330 ± 0.014 (46[. 



This gives xi(2GeV) = (343.8 ±51.6) x 10" 
Xo(2 GeV) = (253 ± 68) x 10~ 6 . 



B. Low-energy theorems 



' GeV" z and 



,^ 2 , 1 f°° , tlmU Q {t) 



(24) 



Unitarity implies that the spectral functions are positive 
for t > t + and satisfy the inequalities 35, 38]: 

Ml(t) ,^I(^±>|z^ l/+(i)|2 , (25) 



2 48tt 



Imn (i) > - 



3M- [(t-t + )(t-t^/ 2 



2 16tt 



\Mt)\ 2 , (26) 



with t± = (M K ± M^) 2 . 

In the limit Q 2 >> A.q CD , the correlators Xi(Q 2 ) 
and Xo(Q 2 ) can be calculated by perturbative QCD. Re- 
cent calculations to order af (see 0, [45j and references 
therein) give: 



xi(Q 2 



a 



8ir 2 Q< 



(1 + — - 0.062a 2 



Xo(Q 2 



3(m s - m u ) 2 
8tt 2 Q 2 



0.162a 3 -0.176(4), (27) 



(1 + 1.80a s + 4.65a 2 

+ 15.0a 3 + 57.4a^). (28) 



We omitted the power corrections due to nonzero masses 
and QCD condensates, as they are negligible. 

The relations (|2"Tj) - ([2l)|) show that each form factor sat- 
isfies a relation of the type (Q| , where 



P+(t) 
Po(t) - 



3/2 



1 [(t-t+)(t-t-)] 

" 32tt 2 t 2 (t + Q 2 ) 3 

3M- [(t-t+)(t-t-)] 1/2 

32tt 2 t 2 (t + Q 2 ) 2 



(29) 



1 The choice of the renormalization-group invariant correlators is 
not unique. Alternative definitions and the corresponding disper- 
sion relations are shown to lead to almost identical constraints 
for the Kis form factors |4:J . Here we consider for convenience 
only the correlators i21l and (122 D . 



The theorems based on symmetries at low energies pro- 
vide useful ingredients in the applications of the above 
formalism. At t = by construction one has /o(0) = 
/+(0), since f-{t) is regular at t = 0, and SU(3) sym- 
metry implies /+(0) = 1. Deviations from this limit are 
expected to be small [l6[ and have been calculated in 
chiral perturbation theory [l7l J47|and more recently on 
the lattice (see the reviews in [l(-[a|). In the case of the 
scalar form factor, current algebra relates the value of 
the scalar form factor at the CT point Ak -k = ME- — M 2 
to the ratio Fk / F n of the decay constants [3 [l9| : 

fo(A Kn ) - Fk/F^/ + A CT . (31) 

To one- loop in ChPT in the isospin limit Act = —3.1 x 
10 -3 |2l|. Results on higher-order corrections, and also 



beyond the isospin limit, are available 0, 
At Ak-k(= —Ak-k), a soft-kaon result 
value of the scalar form factor to F v /Fk 

fo(-A K7r ) = F n /F K + Act- 



relates the 



(32) 



A calculation in ChPT to one-loop in the isospin limit [21( 
gives Act = 0.03, but the higher order ChPT corrections 
are expected to be larger in this case. The estimate made 
in [3(3] leads to a rather large allowed interval for Act- 

In the present work we use as input the values of the 
vector and scalar form factor at t = 0. For the scalar form 
factor we impose also the value Jq{Aktt) at the first CT 
point. As discussed in due to the poor knowledge 
of Act, the low-energy theorem (|32|) is not useful for 
further constraining the shape of the form factors at 
low energies. 

For generality, we shall present our results for arbitrary 
values of the parameters /+(0) and /o(A^-„-). For the 
numerical illustration of the results we shall use as default 
the values 

/+(0) = 0.962 ±0.005, f Q {A Kn ) = 1.193 ± 0.006. (33) 

The central value of /+(0) coincides practically with the 
ChPT prediction given in [TtJ and is quoted in [l| as the 
most recent lattice result. A slightly different average of 
lattice results /+(0) = 0.959 ± 0.005 is also quoted in 
Q. The value of / (A & ) is consistent with the values 
reported in the review I'll in the isospin limit. 
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The value (f5T|) at the first CT point is quite pre- 
cise in the standard model(SM). On the other hand, it 
has been suggested that deviations from ChPT at both 
CT points would be a signature for physics beyond the 
SM, such as right-handed quark couplings to and 
charged Higgs [2, [H, In what follows, we shall de- 
rive constraints on the expansion coefficients both with, 
and without the condition at the CT point. We shall also 
derive a relation that correlates the values of the scalar 
form factor at t = and at both CT points, which acts 
as an independent constraint for specific models beyond 
SM. 



C. Phase and modulus on the elastic region of the 

cut 

We recall that the first inelastic threshold for the scalar 
form factor is set by the Kr\ state, and for the vector 
form factor by the state K*tt, which suggests we take 
tin = (I GeV) 2 . Strictly speaking, we must consider the 
inelasticity due to Kwtttt at (0.91 GeV) 2 , but its influence 
is considered weak and may be neglected. Moreover, it 
is known that the elastic region extends practically up to 
the Krj threshold, which would justify the choice t in — 
(1.4 GeV) 2 . The analysis performed in [4f| led to the 
conclusion that this choice overconstrains the system, at 
least for the scalar form factor. Therefore, in this work 
we make the conservative choice t\ n — (f GeV) 2 for both 
the scalar and vector form factors. 

Below t in the function 8{t) entering ^ is the phase of 
the S'-wave of I = 1/2 of the elastic Kir scattering for 
the scalar form factor, and the phase of the P-wave of 
/ = 1/2 for the vector form factor. In our calculations 
we used as default below ti n the phases from [H, [5(| for 
the scalar form factor, and from [13, H(| for the vector 
case (the differences between the two phases were taken 
as an estimate of the uncertainty related to this input; 
see [4(1). We recall that, while the standard dispersion 
approaches require a choice of the phase above i; n , the 
present formalism is independent of this ambiguity fill ] . 
Above ii n we have taken d(t) as a smooth function ap- 
proaching 7r at high energies. As we mentioned in Sec. [TTl 
the results are independent of the choice of the phase for 
t > tj n . We have checked numerically this independence 
with high precision. 

To estimate the low-energy integral in (|5J), we used the 
Breit-Wigner parameterizations of |/+(t)| and |/o(i)| in 
terms of the resonances given by the Belle Collaboration 
24] for fitting the rate of r — > Kttu decay. This leads to 
the value 31.4 x 10~ 5 GeV -2 for the vector form factor 
and 60.9 x 10 -6 for the scalar form factor. By combining 
with the values of 7 + o defined in (|3"0|) . we obtain 

I' + = (312 ±69) x lO^GeV" 2 , I' = (192 ±90) x 10~ 6 . 

(34) 

Note that we use the Breit-Wigner parameterizations of 
\f+,o{t)\ obtained in [24| only for estimating the low- en- 
ergy integral appearing in ([8} . The parameterizations are 



not extrapolated outside the resonance region; therefore 
their analytic properties do not play a role in the present 
formalism, where the analyticity of the form factors is ex- 
actly imposed. We note also that the dependence on the 
modulus information below i; n is very mild, as it is used 
only for the computation of the low-energy part of the 
integral (j4]). In fact, the results depend rather weakly 
on the values of I', , and moreover the dependence is 
monotonic: an increase of these quantities leads to more 
conservative bounds. 

In our earlier work, [40l | we had carried out an analysis 
on the uncertainties to be attached to our determination 
associated with the uncertainties of the inputs. This was 
essential as the rather high choice of t- m of (1.4 GeV) 2 . 
Having chosen for this work the conservative choice of 
(1 GeV) 2 , par tly motivated by the considerations of our 
later work [4J| , we have obtained the new constraints that 
comfortably accommodate the prior results, including the 
uncertainties. Therefore, we consider only the results 
based on the central values of the input on the unitarity 
cut. 



D. Outer functions 

The outer function is defined in (TIB")) . Using (fTUJ) and 
, a straightforward calculations leads to 



w+(z) = 



1 



Vl -z 2 



8V27rt 

(1 + z{-Q 2 )f (1 - zz{t + )f' 2 {l -zz{t-)fl 2 
X (1 - z ~z{-Q 2 )Y (1 + z(t+))3/ 2 (l + z(i_)) 3/2 
for the vector form factor, and 
V3(M 2 K - M 2 



(35) 



w (z) 



VI -«(! + «) 



3/2 



(36) 



16V27rf 

(1 + S(-Q 2 )) 2 (l -zi(t+)fl 2 (1 -zz(L))'/ 2 

X (l-z5(-g 2 )) 2 (l + 5(< + ))V2(l + 5(<_))l/2 ' 

for the scalar form factor. Here, z is the current variable 
and z{t) is the function defined in (flX)|) . The outer func- 
tions for the standard version [35l . 



37, 38] are obtained 



from the above by replacing t- m with t + . 



IV. PARAMETERIZATIONS OF K a FORM 
FACTORS 

The first parameterizations used simple pole models 
describing the ^-dependence of f+(t) and fo(t) in terms 
of the lightest vector and scalar resonances with masses 
M v and M s , respectively: 



A(t) = A(0)^^. (37) 



The increase of the precision of the experiments re- 
quired more effective parameterizations. 
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A. Taylor expansions 

The simplest expressions, adopted in practically all the 
experimental analyses [f|-[l3j], are based on the Taylor 
expansion around the point t = 0: 



K-l 



K-l 



f + (t) = l+J2ck,+t k , /o(«) = 1 <%,„** (38) 



fc=i 



fe=i 



where f+,o(t) — /+,o (£)//+ (0). The first coefficients are 
often expressed in terms of dimensionless slope and cur- 
vature: 



A' 



A" 



Cl 



M 2 



C2 



2 M 4 , 



(39) 



which are also related by A' = (r^ K /6) and A" = 
2M^ + c to the radius squared (r^ K ) and curvature c used 
in some papers [13, H3, [H| . 

The Taylor expansions converge in the disc \t\ < t + 
limited by the first unitarity branch point. Therefore, in 
the semileptonic range Mf < t < i_, the convergence 
is expected to be rather good, with the asymptotic rate 
t-/t + = 0.31. Of course, the convergence becomes poor 
if the expansions are used outside the region. 

At the present experimental accuracy, only a few co- 
efficients Cfc can be determined from the data (common 
choices are K — 2 or K — 3). A theoretical correlation 
between the coefficients would be helpful in fitting the 
data with more parameters. The formalism discussed 
here is a useful tool in this sense. In the next section, 
we shall derive strong correlations between the slope and 
curvature defined in (|39p . Also, we will show how to 
obtain a bound on a suitably defined truncation error, 
which reflects the influence of the neglected higher order 
terms in the expansions (|38[) . 



B. Dispersive parameterization 

Recently, NA48 [H[, KLOE [H, and KTeV QJ] Col- 
laborations reanalyzed their data with a dispersive rep- 
resentation of the Omnes type, proposed in [29|, |3(| . The 
parameterizations of the two Kta form factors read: 



/+(*) = /+(0)exp 
t 



(A+ + H{t)) 



/ (i)=/+(0)exp 



M\ - Ml 



t 

Ml 

(ln[C] - G(t)) 



, (40) 



where A+ is the slope of the vector form factor, ln[C] = 
ln[/o(M|- — M%)} is the logarithm of the scalar form factor 
at the CT point and the functions H(t) and G(t) are 
dispersive integrals upon the phases 5+ o(t) of the form 
factors. 

The advantage of this type of parameterization is that 
it includes information on the analytic properties of the 
form factors in the complex plane and on the phase, 



known (modulo ±7r) in the elastic region t < t ln from 
low-energy Kit phases. However, at larger energies the 
phase is not known, and this introduces an ambiguity in 
the representation. At infinity, the phase approaches a 
constant, whose value depends on the number of zeros of 
the form factor [3(J [5j , so as to ensure the asymptotic 
decrease like 1/t required by perturbative QCD |52j |. 

Actually, the dispersion relations (140[) require only one 
subtraction. In order to display the free parameters A+ 
and C, an additional subtraction was performed, at t = 
and t = A/^tt, respectively. This reduces the dependence 
on the unknown 5(t) above U n , but at the same time 
spoils the asymptotic behavior of the form factors (see 
^53j for a discussion in a similar context). The correct 
behavior can be restored only by imposing additional sum 
rules, which in general are not easy to implement in the 
fitting procedure. 

As noted above, the representations (|4"0l assume that 
the form factors do not have zeros in the complex plane. 
The influence of possible zeros, analyzed in [30j . depends 
on their position in the complex plane. Information 
about the presence of the zeros at low energies is therefore 
important for the dispersive representations mentioned 
here. In Sec. I VIII we shall derive rigorous domains where 
zero values of the form factors are excluded. 



C. z-parameterizations 

A class of parameterizations used alternatively for vari- 
ous weak form factors are based on expansions in powers 
of a variable that maps conformally the t-plane onto a 
disk. For the Kir form factors the method was discussed 
in [38( and more recently was used by the KTeV Collab- 
oration to reanalyze their data [11J] . 

The expansion is actually based on the method of uni- 
tarity bounds discussed in the present paper. Consider 
the standard version, based only on the inequality 
without information on the phase and modulus on the 
unitarity cut. In this discussed at the end of 

Sec. [HI one should replace t- m by t + when I' = I and the 
functions 0{t) and oj(z) defined in © and p^|) . respec- 
tively, are equal to unity. Then, referring for illustration 
to the vector form factor, from (JT2J) one can write the 
representation 



/+(*) = TT ^9kZ k 7 

w + (z) 



(41) 



fc=0 



where 



VT^(l + z(-Q*))3(l-zz(t-))^ 
+ { > 32v^T(l-zz(-Q 2 ))3(l + z(t_))3/ 2 - [ ' 

Here z = z(t), where 



m = 



(43) 
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An advantage of the z-expansion is that it allows one to 
derive a bound on the truncation error, describing the ef- 
fect of the neglected higher order terms in the expansion 
38]. On the other hand, from (g^J) it follows that the 
outer function vanishes at z — ±1, points that by (|4"3")l 
correspond to t+ and infinity in the i-plane, respectively. 
The zeros in the denominator are not compensated auto- 
matically if the sum in the numerator of (|4ip is truncated 
at a finite order. Therefore, the representation (HIT) has 
unphysical singularities at the threshold t + and at infin- 
ity. These deficiencies of the standard z-expansion were 
discussed in the similar case of the Bn form factor in [53[ , 
where alternative z-expansions free of such singularities 
were investigated. Such parameterizations are useful and 
deserve further study also for the Kg^ form factors. 



In Fig. [21 the constraint (|4"4"1) is represented together 
with experimental points from 0, lalEll2-ll5l, where we 
have extracted the corresponding curvature from the con- 
strained fit given in We note that, except the results 
from NA48 and KLOE, which have curvatures slightly 
larger than the allowed values, the experimental data 
satisfy the constraints. We note also that the theoret- 
ical predictions = (24.9 ± 1.3) x 1CT 3 , A" = (1.6 ± 
0.5) x obtained from ChPT to two loops @, and 

A' + = (26.05±g;|i) x lO- 3 , X'l = (1.29±g$) x 1(T 3 [27|, 
and A^ = (25.49±0.31) x 1CT 3 , A'| = (1.22±0.14) x KT 3 
|3l| | obtained from dispersion relations are consistent 
with the constraint: the expression in the left side of 
(]44fl is negative when evaluated with /+(0) = 0.962 and 
the central values of the sloDe and curvature given above. 



V. CONSTRAINTS ON THE EXPANSION 
COEFFICIENTS 

In this section we consider the most common param- 
eterization of the Km form factors based on the Taylor 
expansions (|38[) . This parameterization does not include 
in an explicit way information on the analytic properties 
of the form factors and their behavior on the unitarity 
cut. However, the matematical method reviewed in Sec. 
[IT] allows one to derive constraints on the expansion coeffi- 
cients, which follow from these properties. One improves 
in this way the quality of the expansion, which includes 
in an implicit way additional theoretical information. 

Using as input the value of I' + given in (IM)) and the 
phase and modulus below ti n = (1 GeV) 2 described in 
Sec. HID we obtain from (fl8|) the following constraint on 
the slope A+ and curvature A" of the vector form factor, 
for an arbitrary /o = /+(0): 

/o[( A +) 2 -0-107A' + A';+2.18xl0- 4 A'|+2.98xl0~ 3 (A' + ) 2 

- 1.49 x 10~ 5 A' + + 4.20 x 10~ 8 ] - 4.67 x 10~ 7 < 0. 

(44) 

The numerical coefficients of this relation depend on the 
phase in the elastic region and the coefficient I', defined 
in (|34l) , which gives the last term in the left hand side of 
the inequality. 

As an illustration, for the input value /+(0) = 0.962 
given in (|33| , the inequality (|44)l is represented as the in- 
terior of the smallest ellipse in the slope-curvature plane 
in Fig. [TJ The alternative input /+(0) = 0.959 adopted 
in [l[ leads to practically the same ellipse. For complete- 
ness, we represent in the same figure the domains ob- 
tained using as input the normalization at t = and the 
standard unitarity bounds, without information on the 
phase and modulus (the largest ellipse), and the domain 
obtained by including in the standard unitarity bounds 
the phase on the elastic region, known from the Fermi- 
Watson theorem (the intermediate ellipse). The small 
ellipse is situated inside the other two, which confirms 
that all the constraints are satisfied by the domain de- 
scribed by the inequality (l4"4"l) . 



-■ norm + phase 




.4 I i I i I i I i I i I i I i I 

-0.2 -0.1 0.1 0.2 0.3 0.4 

FIG. 1: Allowed domain for the slope and curvature of the 
vector form factor, using the normalization /+(0) = 0.962 and 
phase and modulus information up to t- m = (1 GeV) 2 . 



The analogous constraint for the slope and curvature 
of the scalar form factor for an arbitrary normalization 
fo = /+(0) reads: 

/o[( a o) 2 -0-059A[ ) A; ) '-3.58x10- 4 A; ) ' + 9.72x10- 4 (A; ) ) 2 

+ 9.64 x 10~ 6 A^ + 3.67 x 10~ 8 ] - 8.05 x 10~ 8 < 0. 

(45) 

In deriving this constraint, we used the value of I' Q from 
(j34j) and the phase and modulus described in Sec. Mil 
For illustration, for fo = 0.962 we obtain the small el- 
lipse in Fig. [3] where, as in Fig. [TJ the larger ellipses 
are obtained using the standard unitarity bounds. These 
constraints are satisfied by the points in the domain (|45p . 

If we include in addition the constraint at the first CT 
point, fcT = fo(A.Kw)] we obtain the inequality 

fo [(Ao f + 0.25A A ' + 21 .6 x 10 _3 A^ + 15.3 x 10" 3 (A^) 2 

+ 2.68 x 10~ 3 Ao + 1.17 x 10~ 4 ] - 10" 3 / / C t(2.67Ao 

+ 21.53Aq +0.23) + 1.16 x 10" 4 /^. T -3.23 x 10~ 10 < 0. 

(46) 
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FIG. 2: The best constraints for the slope and curvature of 
the vector form factor for /+(0) = 0.962 compared with ex- 
perimental determinations. 
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FIG. 3: Allowed domain for the slope and curvature of 
the scalar form factor, using as input the normalization 
/+(0) = 0.962 and phase and modulus information up to 
i in = (lGeV) 2 . 



For the central values given in ((33)) . this domain is the 
interior of the small ellipse in Fig. 01 The large ellipse in 
this figure represents the allowed domain obtained from 
the standard unitarity bounds with the same input at 
interior points, and the intermediate ellipse is obtained 
from the standard unitarity bounds by imposing also the 
phase below t ln according to the Fermi- Watson theorem. 
The small ellipse is situated inside the other two, which 
shows that the slope and the curvature which satisfy the 
inequality (|4"S)> obey also the standard unitarity bounds. 

The above domains were obtained for the central values 
of the parameters in . As shown in [i(J , the inclusion 
of uncertainties in the inputs has the effect of slightly 
enlarging the allowed domains. In particular, there is a 
straightforward dependence of the shape of the ellipses on 



the input I'. These could be subject to uncertanties both 
in pQCD as well as due to the uncertainties in modelling 
the modulus information. However, the resulting ellipses 
simply shrink or expand if I' is taken to be smaller or 
larger. 

We note that the theoretical prediction of ChPT to two 
loops A = (13.97_?;f ±0.4) x 10~ 3 , X' Q ' = (8.07_£;£) x lfr 4 
reported in [f| is consistent within errors with the con- 
straint (|4"5]l with the default input (1551) : for the central 
value of the slope A given above, the range of A ' allowed 
by (HD is (8.24 x 10~ 4 ,8.42 x lfr 4 ). The same is true 
for the theoretical prediction A = (16.00 ± 1.00) x 10~ 3 , 
X'q = (6.34±0.38) x 10~ 4 obtained in [H from dispersion 
relations. 

As concerns the experimental values, Figs. [3] and 0] 
show that the determinations [l3j (we have extracted the 
curvature from the constrained fit therein) and are 
consistent with the phase and modulus information to- 
gether with the normalization (|33[) . but are outside the 
domain obtained when we impose also the CT theorem, 
with the input value given in f|33|) . 

In Fig. [5j we compare the allowed bands for the slope 
A , corresponding to the end points of the ellipses defined 
by the inequalities (|4"5|) and (T4*6l for the same input as 
above, with the experimental determinations. As noted 
already, the slope predicted by NA48 0] is not consistent 
with the SM input at t = and t = Aktt- This conclusion 
is very stable with respect to the input value of Iq, which 
gives the last term in the l.h.s. of (flBf : it turns out that 
only by increasing this value by a factor of almost 5, the 
allowed ellipse is inflated enough as to include the central 
value of the slope from Q . The relation (HBl implies also 
that, keeping Iq and /+(0) fixed, one must reduce the 
input value at the CT point down to /o(Axir) = 1-138 
in order to have the central value of the slope from Q 
inside the allowed domain. 

We end this section on the shape of the form fac- 
tors with a comment on the low-energy theorem (|32p at 
the second CT point. As discussed in Sec. IIIIB1 due to 
the poor knowledge of the correction Act, this theorem 
does not constrain further the coefficients of the expan- 
sion beyond the domain obtained with the input at t = 
and t = Aktt- On the other hand, using these values as 
input , the formalism presented in Sec. [TT] leads to an 
allowed domain for the value of the scalar form factor 
at t = — Aktt- Namely, from the proper input in the 
determinant (1181) we obtain the inequality 

0.177/ CT + 0.086/ CT + 0.523/ 2 - 0.425/ /cr- 
0.608/o/ct + 0.246/ct/ct - 1-92 x 10~ 4 < 0, (47) 

where fcr = /o( — A^tt). This relation correlates the 
values of the scalar form factor at t = and at the two CT 
points and represents a nontrivial anayticity constraint 
for the predictions beyond the SM. As shown in [40| . 
for the input values (133)) . the inequality ((T7]) leads to a 
narrow range —0.046 < Act < 0.014 for the correction 
defined in (1321) , improving the ChPT prediction reported 




A. Isospin breaking effects 

The bounds given above were obtained in the limit of 
isospin symmetry. We recall that in our convention Mk 
and M n are the masses of the charged mesons. As dis- 
cussed recently jH, the isospin breaking effects should 
be taken into account at the present level of experimen- 
tal and theoretical precision. At low energy, at next-to- 
leading order in the chiral expansion, there are strong 
isospin violations at order (m<2 — m u )p 2 and electromag- 
netic corrections at order e 2 p 2 , while the next-to- next- 
to-leading order terms include corrections up to order 
(m d - rn u ) 2 p A . 

In the formalism considered here, isospin symmetry 
was used in the unitarity relations ((25)) and (|26|) . where 
the matrix elements of the K + ir° and K°tt + pairs enter- 
ing the unitarity sum were related by symmetry. The de- 
scription of the form factors on the unitarity cut in terms 
of resonances suggests that the isospin corrections in this 
region are small (for a detailed discussion see [27]). So, 
the unitarity relations (|2"5)) and ([2TJ]) and, consequently, 
the outer functions defined in Sec. IIII Dl conserve their 
form. We checked also that a change of about 0.02 rad of 
the phase below ti n , suggested to represent isospin effects 
[30| . leaves the results practically unmodified. 

If the symmetry is broken, one must use the physical 
pion and kaon masses in the evaluation of the unitar- 
ity threshold and the phase space factors. To illustrate 
this effect, we recalculated the constraints - d3SJ) using 
the same input on the unitarity cut, with the exception 
of M„, taken to be the mass of the neutral pion. For 
instance, we obtain instead of (1461) the inequality 

fo [(A^') 2 +0.24A^A^ + 21.3xlO-X' + 15 - lxlo_3 ( A o) 2 
+ 2.62 x IO^Aq + 1.14 x 10~ 4 ] - 10" 3 / /ct(2.61Ao 

+ 21.23Ao+0.23) + 1.13x 10" 4 /ct-3-21 x 10~ 10 < 0, 

(48) 

where the definition of the slope and curvature is still 
based on the mass of the charged pion, as in (j38|) . The 
small differences between (14TJ1) and (|4"8")) are due to the fact 
that the position of the CT point and its image in the 
z-plane are slightly changed if the pion mass is modified. 

The inequality ([4"5]) can be used to constrain the slope 
and curvature for the K + ir° form factor using the isospin 
corrections for the same form factor at t = and the CT 
point @,H^|. For instance, by increasing simultaneously 
the values given in §3$), /+(0) by 2% and / (A Kl ) 
by 0.029 [1 EH, we obtain from (gHJ) an ellipse slightly 
shifted towards the right by the amount 8X' ~ 0.0007 
compared to the small ellipse in Fig. 2) The shifted 
ellipse leads to curvatures slightly higher for the same 
slope than those obtained in the isospin limit. 
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FIG. 4: Allowed domain for the slope and curvature of the 
scalar form factor, using the normalization /+(0) = 0.962, the 
value /o(Ajf^) = 1.193, and phase and modulus information 
up to ti n = (1 GeV) 2 . 
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FIG. 5: Allowed bands for the slope of the scalar form fac- 
tor obtained from (|45p large band, and (|46[1 narrow band, 
for /+(0) = 0.962 and /o(Aktt) = 1.193, compared with the 
experimental information. 



VI. BOUNDS ON THE TRUNCATION ERROR 

One might ask how accurate is the representation of 
the Ki3 form factors in the physical region by a few terms 
in the Taylor expansion (|38p . It is also of interest to 
know the extrapolation error if the truncated expansion 
is used beyond the physical region up to, say, the CT 
point t = Aktt- To answer these questions one needs 
an estimate of the higher-order terms neglected in usual 
parameterizations. As discussed in previous works, for 
instance (3(| , [38[ , [53| , the technique of unitarity bounds 
allows one to find a model-independent estimate of the 
theoretical error produced by the truncation of the Taylor 
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expansion. 

We assume that the fit provides values of the coeffi- 
cients Cfc, for k < K — 1, in the truncated expansion (1551) . 
Since the series is convergent in the physical region, one 
expects the influence of the higher terms to be gradually 
smaller. We recall that, asymptotically, the convergence 
rate scales as t/t+, where t+ is the convergence radius, 
and this ratio does not exceed 0.31 in the semileptonic 
region. 

Following, for instance, the suggestion made in [53| . 
one can define the truncation error as the first term ne- 
glected in the expansion: 



\c K t 



K 



(49) 



As further suggested in 53], it is reasonable to increase 
the number of terms until the truncation error becomes 
smaller that the experimental one. 

Without additional information, it is in general im- 
possible to estimate the magnitude of higher terms in 
an expansion, even if it is convergent. Fortunately, in 
the present case, using the fundamental inequality (|18[) 
of Sec. and the relation (fT2")) between the function g 
and the form factor F, it is possible to derive upper and 
lower bounds on the next coefficient ck in terms of the 
coefficients c^, supposed to be known for k < K — 1. 

For illustration, we give below the constraint relating 
the value fo = /+(0), the slope A' and the curvature A", 
to the next coefficient C3. For the vector form factor the 
relation is 



fo \ c l + c 3 (-0.045 + 79.2A' + - 4953.4A'{) 
f 1892.5(A' + ) 2 - 3.39A^_ - 2.07 x 10 5 A' + A" + 134.45A" 
+ 6.24 x 10 6 (A" f + 5.1 x 10" 3 ] - 0.051 < 0, (50) 



while the similar relation for the scalar form factor, with- 
out imposing the CT theorem, reads 

fa \ c l + cs(-0.29 + 8.41A - 3321. 5Aq') 
+ 123.2(A[,) 2 - 0.18\' - 2.04 x 10 4 AqAo' + 445.9AJ,' 

+ 2.87 x 10 6 (A ') 2 + 0.025] - 0.0087 < 0. (51) 

The numerical coefficients in these inequalities depend on 
the particular dispersion relation satisfied by the QCD 
correlators, their perturbative expressions, and the infor- 
mation on the phase and modulus on the elastic part of 
the cut. 

As a numerical example, let us take for the vector form 
factor the central values A+ = 25.09 x 10 -3 and A" = 
1.21 x 10~ 3 H3]. Then, for our choice / (0) = 0.962, 
we obtain from (I5TH) the allowed interval 1.79 GeV~ 6 < 
C3 < 2.25GeV~ 6 , which implies a correction at the end 
of the semileptonic region 5f+(t-) between 3.5 x 10~ 3 
and 4.4 x 10~ 3 . 

For the scalar form factor, taking a point X' = 15 x 
10~ 3 and Aq = 0.69 x 10" 3 , situated inside the small el- 
lipses in Figs. [3] and |H we obtain from (T5"Tj) the range 



1.14 GeW 6 < c 3 < 1.32 GeV~ 6 , which implies a correc- 
tion Sfo(t-) at the end of the semileptonic region between 
2.2 x 10~ 3 and 2.5 x 10~ 3 . On the other hand, evaluated 
at the CT point, the term c^t 3 produces a correction of 
about 0.014 to fo(Aj(„), larger than the error quoted in 
Q33p. So, even for a quadratic parameterization of the 
scalar form factor, the extrapolation to the CT point is 
not precise enough for testing possible deviations from 
the SM. Of course, in practical applications one should 
use in ([50]) and (|5ip the values of the slope and curvature 
obtained from the fits of the data. 

We note that the relations ([SO)) and (15T]) are useful also 
if one uses a cubic parameterization in the experimental 
analysis. In this case, the relations provide theoretical 
constraints for the next coefficient C3 included in the fit. 
The truncation error can then be estimated using a con- 
straint on the next coefficient C4, which is obtained in a 
straightforward way from the general relation (|18[) . 



VII. DOMAINS WHERE ZEROS ARE 
EXCLUDED 

The question of the zeros of the form factors is impor- 
tant from theoretical and practical points of view. The 
dispersive representations of the Omnes type require the 
knowledge of the zeros in the complex plane. The ze- 
ros are important also in ChPT, where their presence is 
required in some cases by symmetry arguments [5ll . l54j . 

A study of the zeros of the pion electromagnetic form 
factor was performed in [551 ] . For the Kir form factors, 
the influence of possible zeros in the context of Omnes 
dispersive representations has been analyzed in [301 ] . The 
absence of zeros is assumed in the recent analysis of 
KTeV data reported in [3. 

The mathematical techniques presented in Sec. Ullcan 
be adapted in a straightforward way to the problem of 
zeros. Let us assume that the form factor F(t) has a 
simple zero on the real axis, F(to) = 0. From the relation 
(TT2"]) it follows that g(zo) — 0, where zo — z(to). We shall 
use this information in the determinant condition (|18[l : 
if the zero is compatible with the remaining information, 
the inequality (|18p can be satisfied. If, on the contrary, 
the inequality is violated, the zero is excluded. It follows 
that we can obtain from (TT8"]) a rigorous condition for the 
domain of points zq (or to) where the zeros are excluded. 
First assume that we use as input only the value of the 
form factor at t = 0. Then from (fT5|) the domain is given 
by 



-go 



-go 

Ht f 

-5(to) 2 



< 0. 



(52) 



Here, I' is defined in ©, go is related to the value F(Q) by 
the relation (TT2"]) , and z(t) is defined in (ITU)) . If we include 
in addition the value of the form factor at some point t\ 
(for instance, t\ = Aktt for the scalar form factor), the 
condition reads: 
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-3o 



,9o 



Hti) 2 



9o 



l-2(ti) 2 
z(ti)z(t ) 

i-i(ti)z(t ) 



-.9o 

z(ti)z(tp) 
L-i(ti)z(t ) 
*(*o) 2 
l-5(toP 



< 0. (53) 



To illustrate the method we use the default input /+(0) = 
0.962 and f {A K7r ) = 1.193. Then, from ([52]) it fol- 
lows that in the case of the vector form factors, sim- 
ple zeros are excluded in the interval —0.31 GcV 2 < to < 
0.23 GeV 2 of the real axis, while for the scalar form factor 
the range with no zeros is -0.91 GeV 2 < t < 0.48GeV 2 . 
If we also impose the low-energy theorem (|3"Tj) , with the 
value of /o(Aa'tt) from ([3"3"| . the condition ([5"3"f implies 
that the scalar form factor cannot have simple zeros in 
the range -1.81 GeV 2 < t < 0.93 GeV 2 . The formalism 
rules out zeros in the physical region of the kaon semilcp- 
tonic decay. 

We have also studied the sensitivity of the variation 
of our inputs. The dependence on the parametrizations 
of the scattering phase shifts; the uncertainties in /+(0) 
and /o(Axtt) are found to be imperceptible. On the other 
hand, the uncertainties on the quantity I' + given in 
lead to the following limits for the region where zeros are 
excluded for the vector form factors: —0.28 GeV 2 < t < 
0.22GeV 2 for the maximum value and -0.36 GeV 2 < t < 
0.26GeV 2 for the minimum value. The corresponding 
limits for the scalar form factor upon inclusion of the 
constraint from the CT points from the uncertainties in 
I' Q , quoted in ([34]), are: -1.60 GeV 2 < t < 0.91GeV 2 and 
-2.26 GeV 2 < t < 0.97GeV 2 respectively. 

The method can be easily extended to the case of com- 
plex zeros. Since the functions are real analytic, i.e. 
they satisfy the relation F(t*) = F*(t), the zeros ap- 
pear in complex conjugate pairs: if F(to) = 0, then also 
F(J,q) = 0. We can implement this condition by formally 
setting in (15"3"|) t\ — > t$ and g(z(ti)) — > 0. So, the complex 
domain where zeros are excluded by the normalization at 
t = is given by 




Re t 



FIG. 6: Domain where zeros of the vector form factor are 
excluded, derived from ()54f) . The small domain is obtained 
without including phase and modulus in the elastic region. 




Re t 



FIG. 7: Domain without zeros for the scalar form factor, de- 
rived from (|54[) for /+(0) = 0.962. The small domain is ob- 
tained without including phase and modulus in the elastic 
region. 



one against variations of the inputs. Here, too, the main 
uncertainty stems from those of I' + and The results 
of these are shown in Figs. IHlandlTUl 
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-,9o 
-.90 



-go 



i-s(*5) 2 i- 

l{t* )z{t ) 



-go 



:(i*)z(t ) 



l-5(t*)z(t ) l-5(to) 2 



< 0. 



(54) 



The determinant can be generalized in a straightforward 
way to include an additional interior point, such as the 
CT point for the scalar form factor. The corresponding 
domains are given in Figs. EHHl On the real axis, the fig- 
ures indicate the points where double zeros are excluded. 
The small domains are obtained without including infor- 
mation on the phase and modulus on the unitarity cut. 
As discussed in Sec. HH this case is obtained formally by 
replacing t m by t + . 



We emphasize that our method is able to give in a rig- 
orous way the regions where zeros are excluded, but says 
nothing about the remaining regions. Thus, we cannot 
answer the question whether the zeros are excluded ev- 
erywhere, or a zero must exist. We recall that we ap- 
plied a formalism that exploits a necessary condition, 
©, which follows from (QJ, the Fermi- Watson theorem, 
and the knowledge of the modulus in the elastic region. 
Therefore, the violation of this condition is sufficient to 
ensure that the zero is not allowed. 



A. Phenomenological consequences 



As in the case of the real zeros, the issue of the stability 
of the exclusion regions of complex zeros is an important 



Our results show that the zeros are excluded in a rather 
large domain at low energies. This provides confidence 
in the semiphenomenological analyses based on Onmes 
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FIG. 8: As in Fig. [7] using in addition the input /o(Aktt) 
1.193. 




Ret 



FIG. 9: Domains where zeros of the vector form factor are 
excluded, derived from (|54|) . using /+(0) = 0.962 and phase 
and modulus information up to t\ n — lGeV 2 , when I' + is 
varied within the errors. The small (large) domain is obtained 
with the maximum (minimum) value of I' + given in (|34[) . 



representations, like those proposed in [2jj,l2fJ], which as- 
sume that the zeros are absent. Indeed, a zero of the 
scalar form factor at — O.lGeV 2 , which would distort 
the shape at low energies [30j], is ruled out by our re- 
sults. On the other hand, as shown in [3(|, a real zero at 
— 1 GeV 2 leads already to a phenomenological form fac- 
tor very similar to one with no zeros. This value is close 
to the extremity of —0.91 GeV 2 of the range without ze- 
ros for the scalar form factor, and the inclusion of the 
additional constraint at the CT point rules out a zero 
even at —1 GeV 2 . In the case of complex zeros as well, 
the allowed zeros are rather remote to produce visible 
effects: for instance, a pair of complex zeros located at 
to = (0.1 ± 2i) GeV 2 , considered in [3(|, is ruled out by 
our results. 

In the case of the vector form factor, the analysis made 
in 30] using data from r decay concludes that complex 
zeros cannot be excluded, due to the lack of information 
on the phase of the form factor in the inelastic region. 
In contrast, as shown in Fig. [5J the formalism presented 
here leads without any assumptions to a rather large do- 
main where complex zeros are excluded. The stability 




Ret 



FIG. 10: Domains where zeros of the scalar form factor are 
excluded, derived using /+(0) = 0.962, /o(Aktt) = 1.193, and 
phase and modulus information up to t ln = lGeV 2 , when 
Jo is varied within the errors. The small (large) domain is 
obtained with the maximum (minimum) value of I' given in 



USD. 



analysis shows that, although for the maximum values of 
I', and Iq from ([34]) the domains do shrink, the implica- 
tions for the phenomenological analyses are unchanged. 



VIII. CONCLUSIONS 

In this work we have studied the shape of the scalar 
and vector form factors in the Kg^ domain, crucial for 
the determination of the modulus of the CKM matrix 
element \V US \. We applied a formalism proposed in [42j . 
which develops the standard method of unitarity bounds 
by including in an optimal way the phase and the modu- 
lus on the elastic part of the unitarity cut. The formalism 
is also suitable for including information from soft-meson 
theorems at points inside the analyticity domain. 

The work reported here is a detailed application of the 
techniques explored in [401 ] , also wiring in the constraint 
to ti n explored in 41| , to the phenomenology of the kaon 
scmilcptonic decay. It uses the powerful modified formal- 
ism for unitarity bounds and constraints from low-energy 
theorems. It goes beyond explorations of [39| which does 
not use phase and modulus information, or j37} which 
does not use modulus information for the scalar form fac- 
tors. Here, the vector form factor is analyzed using the 
modified formalism which has never been done before. 
The modified formalism is employed to isolate those re- 
gions of the real energy line and complex energy plane 
where zeros are forbidden. Thus, this work represents a 
powerful application of the theory of unitarity bounds, 
which relies not so much on experimental information, 
but on theoretical inputs from perturbative QCD. It pro- 
vides a powerful consistency check on determinations of 
shape parameters from phenomenology and experimen- 
tal analyses. It may be noted that in the context of 
B — » D*lv the dispersive bounds obtained in (5|| were 
successfully used by experimental groups studying the 
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decay. 

The method exploits analyticity and unitarity, but dif- 
fers in several respects from the usual dispersion rela- 
tions. In applying these relations, some assumptions 
about the form factor above the inelastic threshold are 
necessary. Moreover, in the Omnes-type representations 
it is assumed that the form factors have no zeros in the 
complex plane. No such assumptions are necessary in 
our approach. Instead, one exploits a dispersion relation 
for a QCD correlator, which is calculated perturbatively 
in the Euclidian region and is related by unitarity to the 
modulus squared of the form factors in the Minkowskian 
region. Positivity of the spectral functions then leads 
to an integral relation of the form (Q} for the modulus 
squared of the form factor, which is the basic relation 
of the approach. From such a relation one cannot make 
definite predictions for the values of the form factors or 
their derivatives, but only derive bounds on these values. 
On the other hand, a remarkable feature is that an arbi- 
trary number of values can be included simultaneously, 
corresponding mathematically to the so-called Meiman 
problem. Thus, the formalism is very useful for finding 
correlations between the values of the form factors at 
different points and for testing the consistency of inputs 
known from different sources on different regions of the 
complex plane. 

In this work, we have focused on the phenomenological 
consequences of the formalism. We have considered the 
standard parameterization of the scalar and vector form 
factors in the physical region, based on the Tay- 
lor expansion at t — 0, and derived constraints on the 
coefficients of the expansion. The results for the slope 
and curvature are given in simple analytic form in Eqs. 
(P ]) -([lg f for arbitrary values of /+(0) and fo(A Kjr ). The 
numerical coefficients in these inequalities depend on the 
dispersion relation satisfied by the QCD correlators, their 
perturbative expressions, and the input phase and mod- 
ulus on the elastic part of the cut. The sensitivity of the 
coefficients to the uncertainty of the input is quite low, 
as shown in [40l |. 

The constraints (|4"4")) - (r4*6"| can be used in the experi- 
mental fits with quadratic polynomials, or for testing a 
posteriori the consistency of the fitted parameters with 
theoretical information from regions outside the 
range. For illustration, the small ellipses in Figs. [TJ3] 
represent these domains for the default input ([33]) . The 
allowed values also satisfy the standard bounds and the 
phase condition in the elastic region, being included in 
the larger ellipses in the figures. 

A more general condition, which correlates the coeffi- 
cients of a cubic Taylor expansion, is given in Eqs. ((50)) 
and ([ST]) for the vector and scalar form factors, respec- 
tively. These relations can be used either to estimate 
the truncation error of the quadratic expansions, or for 
constraining the fits based on a cubic parameterization. 

We have worked in the isospin limit, but have also 
given a brief discussion of isospin breaking effects. The 
relation (|48l) , obtained using the mass of the neutral pion 



instead of the charged one, allows one to correlate the 
isospin corrections in the slope and curvature of the scalar 
K + ir° form factor to those at t = and t = Aktt- 

We have also considered two other applications of the 
formalism: in (|47p we have derived a relation between the 
values of the scalar form factor at the first and the second 
CT points, for an arbitrary /+(0). This represents a non- 
trivial analyticity constraint for the predictions beyond 
the SM suggested in HUES!]. We have studied also the 
possible zeros of the form factors, and have shown that 
they are excluded in a rather large domain at low ener- 
gies for both the vector and the scalar form factors. The 
results support the recent dispersive representations of 
the Omnes type, which assume that the zeros are absent. 

Finally, we point out that alternative expansions, for 
instance, in powers of the variable z defined in Sec. IIV CI 
may be more convenient from the point of view of conver- 
gence than the standard Taylor parameterizations (|38l) . 
Constraints on the coefficients of such expansions can be 
derived using similar techniques, and will be investigated 
in a future work. 
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